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Abstract 

We calculate the polarization tensor of charged gluons in a Abelian 
homogeneous magnetic background field at finite temperature in one loop 
order Lorentz background field gauge in full generality. Thereby we first 
determine the ten independent tensor structures. For the calculation of 
the corresponding form factors we use the Schwinger representation and 
represent form factors as double parametric integrals and a sum resulting 
from the Matsubara formalism used. The integrands are given explicitly 
in terms of hyperbolic trigonometric functions. Like in the case of neutral 
gluons, the polarization tensor is not transversal. Out of the tensor struc- 
tures, seven are transversal and three are not. The nontransversal part 
follows explicitly from our calculations. 



1 Introduction 

Investigations of non-Abelian gauge fields at finite temperature which started 
many years ago have been initiated by two fundamental phenomena, the elec- 
troweak phase transition and the deconfinement phase transition happening at 
high temperature. Numerous methods of calculations were developed and applied 
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and important results on the phase transitions and properties of the fermion and 
gauge boson plasma have been obtained. The key element of most calculation 
schemes is the polarization tensor. It gives information about the spectra of 
fields and it is necessary to calculate various thermodynamic potentials beyond 
tree level. This object was investigated in both, continuum quantum field theory 
and lattice Monte Carlo calculations. In continuum theory, in particular, there is 
the necessity of resummations of perturbation series because of infrared singular- 
ities appearing in higher orders. That was implemented in different resummation 
schemes such as hard thermal loop resummations pQ, [21 E], nonperturbative so- 
lution of the Schwinger- Dyson equations [3] -[10], 2PI effective action [HI [12], etc. 
A wide bibliography on these problems is adduced in the review [T3] . 

In the past decade it was realized in both, continuum field theory [211 1221 [231 
[23] and Monte Carlo computations on a lattice [23 [2H] that at high temperature 
magnetic fields of the order gB ~ g 4 T 2 , where g is the gauge coupling constant, 
B is the field strength T is the temperature, are spontaneously created. This 
phenomenon is important for both high temperature QCD and the early universe 
where strong magnetic fields of different kinds had to be present [26]. Hence, 
a comprehensive investigation of non-Abelian gauge fields in a magnetic back- 
ground at high temperature is in order. Although such investigations also started 
many years ago, they mainly dealt with the one-loop calculations of the effective 
potential [H]-[IH] with the goal to study the influence of temperature on the 
instability of the background field related to the tachyonic mode (n = —1) which 
is present in the charged gluon's spectrum = p\ + (2n + l)gB, n = —1, 0, 1, 
where is the momentum component along the field direction. At finite tem- 
perature, the situation with the instability is different because a temperature 
and field dependent gluon magnetic mass can be generated and possibly remove 
the instability from the spectrum. At present, this problem is not finally settled 
although in some approximation it has been considered in Refs. [T9l [221 123] . 
In fact, only few attention was spent to study the radiation spectra of gluons 
or iy-bosons. Frequently, only fragments were considered like the projection to 
the neutral sector or to tree level states or with special gauge choices [22], [33J. 
In addition, the question about the transversality was not finally settled what 
caused a number of problems over the years (see a recent discussion in [20]). 

In our previous papers [23 [30] we investigated the pure one-loop gluon po- 
larization tensor in a homogeneous magnetic background field, first at zero and 
later at finite temperature. In [29J we calculated the polarization tensor for both, 
the color neutral and the color charged gluons, in covariant gauge at zero temper- 
ature. First we determined all possible tensor structure. It turned out that the 
polarization tensor in the magnetic field is not transversal (although the 'weak' 
transversality, i.e., if multiplied on both side with the momentum, holds). The 
corresponding form factors were represented within the Schwinger formalism as 
double parametric integrals with explicitly known integrands. The nontransversal 
structure was confirmed by an investigation of the corresponding Slavnov- Taylor 
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identities [30] . In Ref . [31] the neutral polarization tensor was calculated at finite 
temperature. The additional tensor structures appearing due to the temperature 
were determined and all form factors were represented again in terms of the dou- 
ble parametric integrals and an additional sum because of the temperature. In 
the present paper we make the final step in this series and calculate the charged 
polarization tensor at finite temperature. 

We are faced with two problems. First, we need to find all tensor structures 
Tp V with the property that only p^T^ u p v = holds where p^ is a non-commuting 
momentum because of the magnetic field (for details see below), which to our 
knowledge had not yet been written down in a complete form for nonzero tem- 
perature. The second problem is to represent the form factors in a way expressing 
the weak transversality in an explicit form. For this it was necessary to find in 
the parametric integrals the necessary structures which allowed to integrate by 
parts. This is done and the remarkable property observed earlier that the surface 
terms from the integration by parts cancel just the contributions from the tadpole 
graphs up to the Debye mass term is confirmed in this case too. 

In the following we follow as close as possible our previous paper [31J. In 
section 2 the basic notations are introduced and the tensor structures of the 
polarization tensor are described. In section 3 the calculation of the form factors 
is actually carried out. For instance, in section 3.5 the Debye mass of charged 
gluons in the magnetic field is obtained. The Conclusion summarize the main 
results, the explicit expressions for the one-loop form factors and some discussion. 
The procedure of restoring the form factors from the expressions calculated in 
the main text is given in the Appendix. 

Throughout the paper we put h = c = k B = g = l and also the magnetic field 
B is put equal to one in the most part. The dependence on B can be restored 
by T — > TIB 1 ! 2 f or the temperature and p^ — * p^j ' B 1 / 2 for the momenta. 



2 Basic notations 

We start from the operator structures T^, which are allowed by the weak transver- 
sality condition, 

Px tS p x , = , (1) 

which follows from the corresponding relation the polarization tensor has to obey. 
Here p\ is the momentum depending on the external field. These structures 
appear in the expansion in terms of form factors which will be given below. In 
the magnetic background field, the polarization tensor can be constructed out of 
the vectors / M , and d^, 
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where the third vector is = F^ u p u and we note p x = l x + h\, and the matrixes 
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Hence the operator structures can also be constructed out of these quantities 
only. We mention that with our choice B — 1, in ([3]) is the field strength of 
the background field. 

Further, the structures T xx , can be at most quadratic in the momenta and 
following the polarization tensor they must be Hermite. Writing down all allowed 
combinations, a set of linear independent ones is 



T (l) 

AA' 


— l 2 fi\\> — hh' 






T m 

AA' 


= h 2 5i x , + 2iF xx , 


- h x h X i 


= d x d X i + iF XX i 


T (3) 
AA' 


= h^ + FSty 


- l\h x > - 


h x l x > 


T (4) 
AA' 


= i(l x dy - d x l\') 


+ il 2 F xx , 


-4' 


AA' 


= h 2 5l,-l 2 5i x , 






rT-i(6a) 
AA' 


= 5 l l x ,+l 2 lF X y, 


AA' 


= 3<&, + h 2 iF xx 



(4) 

where also the identity i(d x hy — h x d X i) = ih 2 F xx > +8 XX , was used. In our previous 

"" and T xx ) only their 



paper [29] we mentioned instead of the two structures T AA 
sum, 

rp{6) _ rp(6a) rpi&b) /r\ 

1 AA' — J AA' + 1 AA' • l°J 

In the end it will turn out that for the considered one-loop contribution this is 
sufficient. However, for the calculations in this section it is convenient to keep 
temporarily separately both, Tjffi and T AA ?. 

Another remark on the properties of the structures T< 



AA' 



is that the first four 



are transversal, PxT xx , = T xx ,p X i = holds for i = 1,2,3,4 in addition to ([T]). 
The first three structures are just a decomposition of the kernel of the quadratic 
part of the action, Eq. (24) in [29J, 
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In the case of finite temperature, in addition to the magnetic field, we have 
to account also for the vector u x which describes the speed of the heat bath. In 
the following we assume it to be orthogonal to h x . In fact we use u x = (0, 0, 0, 1). 
With this vector additional tensor structures obeying the weak transversality 
condition ([[]) can be constructed, 



T 



XX' 



(u X ly + l X Uy) (up) - b\y[up) 2 - U X Uyf 
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T \x = { u \h\' + h\U\i) (up) - 5 xx ,(up) 2 - uxuyh 2 , 
T l\> = u x id X ' - id x u X ' + 2iF X y(up) , 

TR? = 5 l l xl (up) 2 -u x uyl 2 , Tj$ = 5i x ,{upf - uxuyh 2 . (7) 

These structures are linear independent. Below, however, it will turn out that 
rpioa an( j rpiob a pp ear a ^ fa e one-loop level which we consider here only in inter- 
mediate steps and drop out from the final result. 
In addition to there exists a further structure, 

= u x u y (8) 

which fulfills (0Q) for p^ = (in an obvious, trivial way). Having in mind that the 
condition ([1]) holds only if the external moment p of the polarization tensor fits 
into the Matsubara formalism, i.e., if it is given by p^ = 2nTl (I integer) than 
(jSJ) makes sense as a structure being present for I = only (or, formally, being 
proportional to <5;,o)- It is just this structure which delivers the Debye mass term. 
It must be mentioned that T xx , is contained as special case also in T xx ,, T xx ,, T^f 
or in T xx , but due to its exceptional role it makes sense to keep it as separate 
contribution. 

Before writing down the expansion of the polarization tensor in terms of form 
factors we mention one property which follows directly from the basic commutator 
relation the momentum p^ obeys, namely pxp 2 = (p 2 8xx + 2iF\y)py. As a 
consequence, for a function of p 2 the relations 

Pxf(p 2 ) = f{p 2 + 2iF)xx>Px> , 

f(p 2 )Px = pyf(p 2 + 2iF)yx (9) 

hold where now / must be viewed as a function of a matrix so that it itself 
becomes a matrix carrying the indices A and A'. The same is true with h 2 in 
place of p 2 . 

The decomposition of the polarization tensor can be written in the form 

U X y(p) = E nW (^ h2 + 2 ^)aa» t§ x , + n D r A D A ,. (10) 

i 

The sum includes in general all structures T xx , defined in 03J) and in . The form 
factors U^(l 2 ,h 2 + 2iF) X y depend on I 2 and h 2 only (besides their dependence 
on the matrices in (jSJ)). In ffTOl the form factors can be placed also on the right 
from the operator structures applying both relations (jUJ). 

3 Calculation of the polarization tensor 

The basic Feynman graph for the polarization tensor is shown in Fig. 5 in [2H] 
and the notations of the vertex factors in Fig. 2 there. The analytic expression in 
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momentum space reads 

/elk 
{T\ up G uu >(p- k)T X 'u'p'G pp '(k) 

+(p - k) x G(p - k)k x ,G(k) + k x G( P -k)(p- k) x ,G(k)} 

+nS po1 , (ii) 

where the second line results from the ghost contribution and the tadpole contri- 
bution is given by 

nl a A ' po1 = Jj^{hyG pp (k)-G(k) xx ,} 

+ / (^4 + G ^(P) ~ 2<W)} • (12) 

The vertex factor is given by 

T Xvp = {k- 2p) p 5 Xv + 5 pv (p - 2k) x + 5 pX (p + k) v . (13) 
For the following we rearrange it in the form 

T Xvp = (p - 2k) \ 5 vp + 2 {pJxp ~ P P hu) + (-(p - k)Jx p + k p 5 Xv ) 



F (l) _j_ F (2) , p(3) 

— A 1 Ai/p 1 Ai/p> 

where T x J p will be temporary further subdivided into two parts, 



(14) 



rS$ = -(p - *) j Xp and rg = v>- (15) 

The momentum integration in the polarization tensor is carried out using the 
formalism introduced by Schwinger, [32]. There the propagators (in Feynman 
gauge, £ = 1) are represented as parametric integrals, 

p OO P oo 

G{p -k)= ds e~ s{p ~ k)2 , G(k) = / dt e~ tk2 (16) 
Jo Jo 

for the scalar lines and 

POO 

G vv ,{p-k) = / dse- sip ~ k)2 E uv , 
Jo 

Gppik) = / dte- tk 5 pp/ (17) 



for the vector lines with 

E w , = e~ 2isF = 5{ l x , - iF xx , sinh(2s) + S^ x , cosh(2s) 
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In this formalism the momentum integration is written as an averaging procedure 
in some auxiliary space and for the basic exponential 



6 = e - s{p ~ k) e' tk (19) 

it holds 



rlk 



(27T) 
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with 



Q{l 2 X)= ex P[~^ (21) 
(4vr) 2 (s + t)\/A 

which is the result of the corresponding calculations, see [32] for details. The 
following notations are used, 

H = —I 2 + m (s,t)h 2 (22) 

s ~\~ t 



and 



m(s,t) = s — arctanh- 

Q 



as well as 



lb 1+ , 2t -\" 2 (23) 
2 1 - (1 - 2t)e- 2s v ; 



(g 2 -p 2 ) I A 

t 2 + tsinh(2s) + p/2 (24) 



with the notations 



p = cosh(2s) - 1 , (25) 
q = 2t + sinh(2s) , (26) 

which will be met frequently in the following. With these notations the self energy 
graph for scalar lines becomes represented by the parametric integrals in the form 

n^ifo) = Jj^ g (p- k ) G ( k ) = J dsdt ^> h 2 ). (27) 

These formulas are derived for T = 0. To include nonzero temperature, within 
the Matsubara formalism we are using we have to substitute the integration over 
the continuous momentum by a discrete sum over I in k& — 2irlT, 

P - T E • (28) 

1= — CO 
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In order to incorporate this into into the parametric integral we represent 

oo oo roo 

T = T £ dk 4 5(k 4 - 2%IT) 

i I J —oo 



l=—oo l=—oo 
oo 



00 /■oo roo AU „ 

T Y da e-^ lT / ^ e^ 4 (29) 

7 j —00 J —00 2tt 



in this way keeping the original formalism on the expense of accommodating 
the additional factor exp(iak^) and carrying out the integration over a and the 
summation over I afterwards. Under the above assumption that the speed of the 
heat bath is orthogonal to h\, the additional factor exp(iak4) can be incorporated 
into Schwinger's formalism quite trivially because the integration over k^ (in the 
same way as that over k^) decouples from the other ones in the sense that the 
corresponding integrals factorize. For the integration over fc 4 we have 



OO 



2tt 



and 



T V [°° da [°° e"*** e-d*-**) 3 -^ 
L ' -00 J-oo 2tt 



l=—oo 



00 / 2 y- 

= T V exp -(2vrZT) 2 + 47r/TsM - —(up) 2 - — — 

Z= co V S + * S + * 

In this formula we used (up) = P4 because this form of writing will be more 
convenient below. 

Combining these with (120]) . (12TT) and (120], we get for a graph consisting of 
scalar lines 

00 . j3l 

n scalar s T £ J j—^ G(p-k)G(k) (32) 

OO r 2 

= TY, j dsdt e -( 2 ^) 2 ( s +*)+ 4 ^M-feW 2 ^ 47r(s + t) ( Z 2 j ^ 
i=— 00 

This representation is still not in a form which is sufficiently convenient for the 
following, for instance, it contains still the ultraviolet divergence which appears 
from small s, t and large I. Using the well known resummation formula 



exp C— zl 2 + aPj = J — E exp ( - 

l ' Z N \ 



71 2 N 2 . a a 2 \ 

+ miV- + — (33) 



z z ' Az 2 j 

(both sums run over the integers) we obtain with z — > (2irT) 2 , a — > <AitT(up) 

n scalar = ]T | dsdt 6 T (/ 2 , /i 2 ) . (34) 



N 
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Here we introduced the basic average 



Q T (l 2 ,h 2 )= exp|- 4(g ^ )T2 +2 g (np)|e(/ 2 ^ 2 ) , (35) 

where 0(Z 2 , h 2 ) given by Eq.( l2TT) and where the notation 

iN . . 

ma = —, r^; u\ (36) 

2(s + t)T v ; 

was introduced. This average is what comes at finite temperature in place of 

(EQD, 

T £J$F =(Q)T = QAi 2 ,h 2 ). (37) 

In (I35p the ultraviolet divergence is in the N = contribution taken at B = 0. 
The well known basic properties of the representations as sum over I or as sum 
over N are that the sum over I is convenient for high temperature (/ = - the 
so called static mode, gives the leading contribution for T — > oo) and the sum 
over N gives the low temperature expansion, for instance the N = term is the 
contribution at T = 0. 

At this place it is meaningful to show how the magnetic field can be restored 
in the expressions involving the parameters s and t. Since in the proper time 
representation of the Green's functions fflBT) . ffTTj) the phase is dimensionless and 
since the restoration of the 5-dependence has to respect that, we must take 
s(p — k) 2 — > (sB)((p — k) 2 )/B. Hence in the parameters s and t the magnetic 
field is restored by s — > sB and t — > tB. 

We continue with the remark that since below Q(l 2 ,h 2 ) will become part of 
the form factors it is meaningful to write it as a function of h 2 + 2iF. This can 
be done by means of the relation 



with 



e(l 2 ,h 2 ) =Q(l 2 ,h 2 + 2iF)Z (36 



2 - ~ ET w - " + w iF + w s± • < 39 > 



where the notations 



and 

are introduced. 



a = [p 2 + q 2 ^j sinh(2s) — 2pgcosh(2s), 

p = (p 2 + g 2 ) cosh(2s) - 2pgsinh(2s), (40) 

A = E - 1 , D = A — 2itF (41) 
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From Schwinger's formalism we need also the commutator relation 



Pfi e(l 2 ,h 2 )=e(l 2 ,h 2 % (42) 

with 

V, = (Ep)^ - (Ak)^ . (43) 

Here we used obvious short notations like (Ep)^ = E^p^i. Finally, we need the 
average formulas for vectors, 

(6A; A > T = (<d) T (^ p + u 



(Qk x k x ,) T = (8) T 



D 

A \ (A \ (%F 



D P + U )Ad P + U ) x , + \DT, vv 



(44) 



together with the explicit representation 



± = -?-SW- -^iF + y + tsinh(2s) S\ (45) 
D s + t 2N 2N K J 

All other quantities can be calculated from these, for example, 

iF 1 /-2iF\ T 1 / A\ T 1 / 5 11 p ._ q _,\ 
W = 2{—) =2t( l -D) = 2\—t + jN tF+ jN S ) ■ (46) 

Perhaps it is useful to remark that all these matrices commute one with the other, 
that the transposition changes the sign of F and that the simple algebra 5" =5", 
5 1 - 2 = 5 X , F 2 = -6\ 5^ = S^F = and 5 ± F = F holds. 

The averages in (1441) are calculated at T ^ 0. For T = they reduce to the 
formulas known from [32]. For T ^ one needs to consider the corresponding 
generalizations of Eq. (l3~0]) . 

oo Ah, n „ 

^ fc 4 e iCTfc4 e -(w-**) a -t*J . (47) 

-oo 2lX 

Replacing the additional factor fc 4 by i-^, after integration over a an additional 
factor 27r/T appears in a formula which generalized Eq. (|32|) . It remains to do the 
resummation from I to N. Taking the derivative with respect to a from Eq. (l33"l) 
after some calculations the first line in (144)) for A = 4 appears. The derivation of 
(144")) is then finished by the remark that for A = 1, 2, 3 no additional contributions 
appear. In a similar way also the second line in Eq. (144 j) can be derived. 

Now we turn to the calculation of the polarization tensor (1111) . Using (j!7p . 
flUD, (EHD and it can be written in the form 

fI A A' = J2 J dsdt (® r ^\ypEvv'^\'v'p + (p - k) x k x > + k x (p - k)\>]) 



N 



+ n A a A 1 po1 , (Ah 
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where in Y\ vp one needs to substitute p by p. 

At the next step we divide the whole expression into parts according to the 
division made in (fT4l) and (fl5i) . 



Hav = E / M ( r 

AT ^ 



)+ n 



tadpol 
AA' 



with 



r(>) p pO') 
1 Xvp-^w' 1 X'v'p 



(49) 



(50) 



and Mfx' * s the corresponding contribution from the ghost loop. 

The sums in (j4"9~|) include also the decomposition (|T5|) . The explicit expressions 
for these quantities read 



M 1 ' 1 
M 1 ' 2 
M 1,Sl 

M 2 ' 2 

M 2 ' 31 
M 2 ' 32 

M3 1)2 

m 32 ' 2 



{p-2k) x {p-2k) x , tiE, 
2(p-2k) x ((E-E T )p 



-(p-2k) x (E (p-k)) x ,, 
Z(p-2k) x (E T k) x/ , 
2((E T -E)p)^p-2k) y , 



45 X \> (pEp) 



4 (E T p 
4 



+4(Z) xx ,(pp)-4(p) x , (Ep) x , 
-25xx> (pE(p-k j) + 2p y (E(p 
-2{E) XX , (pk) + 2 (E T p) a# fc A , 

-(£ T (p-A;)) A (p-2A : )v, 



-*)) 



A ' 



(Ek) x (p-2k)y, 

-26xx> {{p-k) Ep) + 2 (P T (p - fc)) ^ p A , 
-2 [E) xx , (kp)+2k y (Ep) x , 

6xx> ((P-k) E(p - fc)) + P AV P - A; a (E T (p - k)) y - (E (p - fe)) A A; A 

(j)-k) x k x > + k x (p - 

(51) 

Using Eq. ( )43|) and simple relations like A + 2 = + 1 and EE T = 1 these 
expressions can be simplified further. After that we apply the average formulas 
QUI and pass from 6(/ 2 , h 2 ) to 0(/ 2 , h 2 + 2zF) by means of (1381) which brings a 
factor Z to the M iJ which couples to the index A, in detail, ZM !J stands for 
Zxx"M l x i x ,. In this way we come to 



(6M !J ) T = 6 T (/ 2 , h 2 + 2iF) M 



1, j 



(52) 



with 



M 



1,1 



(Pp - 2m) A (P T p - 2u 



2Z 



2iP 

p^ 



trP 
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= (Pp-2u) x (Rp) x ,, 
= (Pp-2u) x (-Q 31 p + u) x ,-ZE 



M 1,2 

M l,3i 
M l,3 2 

M 2 ' 1 = (i? T p) (P T p-2u) / 



M 2,2 

M 2 ' 31 
M 2 ' 32 

M 3i,l 
M 3 2 ,l 
M 3i,2 

M 3 2 ,2 

M 3 ' 3 + M gh 



2iF 



(Pp-2u) x (Q 32 p + u) x ,-ZE 



2iF 



-4 



(5p) A (Tp) A , + (T T p) x (s T p) y - S\\i ( P T T p) - T xx , (pSp) 



+AST T 2iF , 



= 2 (V T P - u) x (S T p) x/ - 2Z XX> {(p-§rp) - (<ip)) - 2iFZ 

= 2{Up + u) x (Tp) x , - 2 (EZ) XX , ( (i) T + + 

= {-QliP + u) x (P T p - 2u) a , - Z£ 

= (QI2P + u) x (P T p - 2u) y - ZE 



2%tF 



2iF 



2iF 



(53) 



2iF 
2itF 



= 2 (Sp) x (Vp - u) x , - 2Z XX > ((p-§tp) - (up)) - 2iFZ 



,2itF 



2(T T p) x (u T p + u) xi -2(EZ) xx ,^l [ p 
A 



D 



1 - 



D 

+ (EZ)xx> ( (^P + fi) P + <i) + tr^ ) , 



p — u 



l-^)p-u) ) -U:E 



pj + (-up) 



ZA 



2iF 



where in the last line some cancellations occurred. Here, again, new notations 
were introduces, namely 

P = (1-2^ T = Z(F-(F + 1)^), 

Lll S \\ - % iF - P/l^f. 5 1 - = r 3 6" + a 3 iF + (3 3 5 ± , 



R = 



s + t " A A 
2(e~ E t ) = -4sinh(2s) iF , 



Qsi = 



= E- 



-2UF 



t ^ ^pcosh(2s) — gsinh(2s) ^ ^gcosh(2s) — psinh(2s) ^ ± 



s + t 2 A 

s 4 5 11 + ^ 4 iF + 646^, 



2A 
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Q 32 = (^) T = ^, 



s ^11 p(sinh(2s) + t) pcosh(2s) + tsinh(2s) 



s + t 2A 2A 

= s 3 5 11 +^ 3 iF + S 3 5 ± , 

S = Z=(E-AiY = 6\\ + ±.iF + JLv 



DJ 4A 4A 

= ri^ 11 + aiiF + PiS ± , 

T = (EZy = l+A^ = 5" + 2 ^iF + £±^6\ 

D 4A 4A 

= s 2 5" + 72 i-F + ^2 

V-D/ s + t 2A 2A 
= ratf" + a 2 iF + /3 2 5 ± , 

^ D S + t " + 2A + 2iV " ' 

= s 1 ^+ ll tF + 5 1 5 ± } (54) 

where ( |39|) . (14"5]) and (|46|) were used. The notations rj «i, Sj, 7^ and 5j are 
introduced here for later use. 

As for the dependence on u we made use of the fact that only the fourth 
component of w M is nonzero such that relations like Zu = Su = Tu = u and 
Ru = hold. 

In order to continue and, for instance, to find the necessary structures for 
integration by parts, we divide the contributions into 3 parts, 

1. M 1 ' 1 + M 1,2 + M 2,1 + M 2,2 , 

2. M 1,3 + M 2 ' 3 + M 3,1 + M 3,2 

3. M 3 ' 3 + M gh 

and consider them individually in the following subsections. 

3.1 Contribution from M 1,1 + M 12 + M 21 + M 2 < 2 
We start with M 1 ' 1 . First of all we note that the relation 

-(5-i)(*-^)-»£ < 55 > 
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holds which can be verified by differentiation of (|39|) and (146]) . It should be men- 
tioned that the term — 2 " A -"Y in the left hand side vanishes under differentiation. 

(up) 

It was added by hindsight in order to carry out the integration by parts and to 
derive a convenient representation for the form factor. 



Eq. (155j) allows us to represent the contribution of M ,x to the polarization 
tensor (|49p in the form (up to the sum over N) 

J dait B T (fl, + «F) (- g - |) (P - 2«™ ) ) *E . (56) 

In this integral we temporarily change variables to A and u according to s = Aw, 
t = A(l — u) and with (Jj — ^) = i^: we integrate the derivative with respect 
to u by parts, 

j~ d\\ £ du 6 T (! 2 , h 2 + 2iF) Ll±(p- 2 50) \ tvE 



"rfAA -G^ 2 . h 2 + 2iF)(p - 2^-) 



o 



-e T (l '\ h 2 + 2iF) P - 2-^4 tr£ 
V M J 



u=0 



+ J dsdt 6 T (/ 2 , h 2 + 2zF) j \^P - 2 ^yj ((pPp) - 2(up)) trE 

+ 4sinh(2s)pl , (57) 



where 

<*»~(!-i)*-E(5-s) A <" 



follows from fl22lh (124] and fl54j) and 



(59) 



from (1351) . We have to take into account that we differentiate a G, which depends 
on h 2 + 2iF, so that we get in place of (1581) 

{pPp) = {pPp)\ . (60) 

In such expressions we shall use the tilde as notation for this substitution in the 
following several times. Further we used trE = 2(1 + cosh (2s)). 

As a result from the integration by parts we get surface contributions which 
have the general form 



M surface = / dXX 6 r (/ 2 , h 2 + 2iF) ]T M s 

i,j 



(61) 

u=0 
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The surface contribution from M ' x is then 



1.1 

sf 



P\\> - 2^ I trE . 



(up) 



(62) 



Taking the contribution from the last line in ( 1571) we represent up to the 
surface contribution M 1 ' in the form 



M 1 ' 1 = l^Pp - 2u) x (P T p - 2u) y - ( P AV - ) (ipPp) - 2(uv) ) ! trP 



+4 sinh(2s) P ; 



XX' 



' (up) J 

Next we consider M 1,2 + M 2,1 . It is useful to rewrite these in the form 
M 1 ' 2 + M 2 ' 1 = (Pp - 2m) A (Rp) x , + (P T p) a (> T p - 2«) , 



(63) 



AA' 



' (up) 



(pPJp) - Rj x , ((pPp) - 2(up) 



+ Px.-2^' 



(up) 



-4sinh(2s)) + Rj x , ((pPp) - 2(up))(64) 



In the last line we use (pR T p) = — 4sinh(2s). In fact, by the next-to-last and the 
last lines we added zero. Now we shall integrate by parts in the last term in the 
last line. Using the same procedure as before and, for instance, Eq. (l51?|) . we get 
the surface contribution 



and using 
we represent 

M 1,2 + M 2,l 



M^ = -e T (i 2 ,h 2 + 2iF)R T xx , 

(£-|)l& = 8cosh(2a) zP 

(Pp - 2u) x (Rp) x , + (R T p) x (P T p - 2u 
- (p\\i - 2^f) (pR^p) - Rj x , ((pPp) - 2(up) 



(65) 
(66) 



(up) J 

-4sinh(2s) (P A a' -2^4^- J +8cosh(2s) iF ; 



(up) J 

Finally we turn to M 2,2 . Here we rewrite 

(pSp) = (pSp) + (S + S T ) iF , 
(pT T p) = (pfTp) + (T + T t ) iF 



XX'- 



(67) 



(6£ 
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and represent M 2,2 in the form 

-4 [{Sp) x [Tp) x , + (T T p) x (s T p) y - S xy (pT^p) - T x \,{pSp) 
-8cosh(2s) iF . (69) 



M 2 ' 2 



The last line is the result of a number of cancellations. 

We continue with the observation that the last lines in ( 1631) . ( 1671) and ( l69i) 
compensate each other so that we are left with the corresponding first lines. 
These are written in a form as considered in Appendix A so that we can apply 
directly the equations H108[) , f 1 11 3 j) , (11161) and for the temperature part (11191) and 
so on. The first line in Eq. (l63l) has the form as given by (11141) and (11151) resp. 
(11261) and (\127\\ multiplied by trE, the first lines of Eqs. (1671) and ( 1691) match the 
structure of W 3 in (11091) resp. fll 18j) . The parameters r { a%, Pi, Sj, 7, and 5i are 
given in Eq. (l54|) . 

In this way we get from this subsection the following contribution to the form 
factors, which we denote by Mf: 

Mf = - (f^f) 2 2(1 + cosh(2s)) +8, 

M a = (f P ) 2 -(p/2-f 2 ) 2 2{l + cosh ( 2g )) -8sinh(2s)| +8cosh(2s), 

Ml= f=f^ 2(l + cosh(2s)) +^±2!, ( 7 °) 

M t = -frfl 2 ( 1+ ch ) +4sinh(2s)f=| +^ 

(11) (12) + (21) (22) 

where in the last line the origin of the contribution is indicated. We note that 
there are no contributions to M 5 , M 6a or M eb . 

For the temperature induced contributions we get with Eq. (11271) with fi = 

2 iN 



'2(s+t)T 



^ ^ „, „ iN (3; 



MS = -2—^-— ~r~~r trE, MZ = -2—^-— -Ar trE 



7 2{s + t)T {up) ' 7 ~2(s + t)T {up) 

iN 
'2(s + i)T 



M° = -Mfo = 2 ^- ; ^ a 3 trE. (71) 



3.2 Contribution from M 1 ' 3 + M 2 ' 3 + M 3 ' 1 + M 3 ' 2 

In this subsection we employ the subsplitting (|T5|) . First we consider M 3l,1 +M 31 ' 2 
from (1531) and do a reordering using ([5] 



M 3i,i + M 3 l5 2 = 2 (Sp) A (Fp - u) v -2Sxx' ((pV T p) -up) -2(v + V T )SiF (72) 



\ / \ / \ 2iF 2HF 

-QIP + fi) A ( pT P " 2fi) A , + 2 (V + V T ) SzF - ZE ~y - ZiFZ— 
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The reason for the reordering is that the last three terms in the second line collect 
into a derivative, 

2 (V + VA SiF - ZE^-% - 2iFZ^ = ( #- - - *™) , (73) 

V / D T D T \ds dtj V^ 31 (up) J 

a relation which can be checked by explicit taking the derivatives. For instance, 
the term — vanishes under differentiation. It was added by hindsight. In a 
very similar way we rewrite M ' 3l + M 2,Sl in the form 

M i,3 1+M 2, 3l = 2 (V T p - u) x (S T p) xi -2S X \> ((pV T p) -up) -2(y + V T ) SiF 

+ (Pp - 2fi) A (-Q 31P + fi) A , + (| - |) (flj - ■ (74) 

Now we integrate by parts in the same way as before. The surface contribution 
is 

M s 2 f = 29 r (/ 2 ^ 2 + 2zF) (QJx-^) (75) 
and in the bulk part we get 

M 3x,l + M 3i,2 + M l,3i + M 2,3x 

= 2{ (Sp) A (Vp - u) x , + (y T p - u) x (S T p) y - 2S xy ((pV^p) - up) } 

-{(Pp-2u) x (Q 3l p-u) x , (76) 

+ (Qj lP - u) x (P-P - 2u) x/ - 2 ( Q l - {(pPp) -up)}, 



where also relations like (tool) were used. We observe that the two lines in the 
last formula have just the structure as one of that given in Appendix A. For the 

iN 
2(s+t)T 

iN i ,, _ iN 



upper line it is W 2 with /i = and v = , ^ )T multiplied by an overall factor of 2 



and for the lower two lines it is W\ with \i = ^ 2 {s+t)T anc ^ v = 2{l+t)T multiplied 
by an overall factor of (—1). The corresponding contributions to the form factors 
can be obtained by means of (jIDgj) . (TTT9D and (TT23D . 

In a similar way the contributions from M 32 ' 1 + M 32 ' 2 can be written as 



M 3 2i i + M 3 2 ,2 = 2(T T p) A (f/ T p + M) A/ -2T A T v ((pf/p) + np)-2(f/ + f/ T )T T iF 
+ (QJ2P + u) x (P T p - 2u) y + 2(U + U T ) T T iF - Z^ (77) 
and and that from M 1,32 + M 2 ' 32 are 



M i,3 2+M 2,3 2 = 2(Up + u) x (Tp) x , -2T x T x ,((pUp) +up) ~2(U + U T )T T iF 



+ (Pp - 2u) x (Q 32 p) x , + 2(U + U T ) T T iF - Z^. (7 
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In this case the derivative has the form 



2( U + U ^_^ = (°°)U + ™) (79) 



D T \ds dtj V (up) J 

resulting surface contribution is 



where again is a term not contributing here but added by hindsight. The 



M 2 = -2Q T (l\ + 2zF) [Q] 2 + (80) 
and the corresponding bulk parts after partial integration read 

M 3 2 ,l + M 3 2 ,2 + M l,3 2 + M 2,3 2 

= 2{ (Up + u) x (Tp) x , + (T T p) x (u T p + «) a# - 2T A T V ((pUp) + (up)) } 
+ (Pp - 2u) x (Q 32 p + u) x , + (Qj 2 p + u) (P T p - 2u) 



2 \Ql + ^yj ((PUp) + (up)) . (81) 



These have also the form as given in Appendix A, in both with W±. In the upper 
line one has to take /i = — 2 (s+t)T an< ^ v = ® w hh an overall factor of 2 and 
in the lower lines /x = 2 2 ^+t)T an< ^ v = ~ 2{s+t)T • Their contributions to the 
form factors can be obtained by Eqs. (11081) and HI 131) and fl 1 191) and again we 
use the notations r, «i, Sj, 7, and 5i taken from (J54"j) and r 4 = r 3 , a 4 = ct 3 , 
/?4 = /?3. Further we take into account the numerical prefactors like 2 and -1, and 
we denote the contributions to the form factors originating from this subsection 
by M\. These read 



M b 2 = (-cosh 2 (2s) + (2-4£ 2 )cosh(2s)-12tsinh(2s 

M\ = , 1 — f-2t 3 - 10st 2 - 2 (t 2 + st + 2) cosh(2s)t 
3 2(s + t)AV V J K J 

-8(s + t) sinh(2s)t + 5t-(s+t) cosh 2 (2s) + s 

M * = 2(s + t)A ( 4t(t ~ S) ~ (g + t} + X ) Sinh(2g) 

+ cosh(2s)(4(s - t)i + (s + t) sinh(2s))) , 

Ml = ((cosh(2s) - 1) (2t 2 + cosh(2s) - l)) , 
Ka = ^((-2t 2 + cosh(2s)-l)sinh(2s)), 

Ml b = — ((-2t 2 + cosh(2s) - l) sinh(2s)) . (82) 
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We observe that M| a and Mg b are equal. Also, we remind that there was no 
contribution to M 6 from the preceding subsection and there will be none from the 
next subsection. As a consequence, the operators structures Tjffi and Tffi come 
with the same form factors, hence these contributions collect into the structure 
T$, Eq.©. 

Finally we collect the temperature induced part. Its form factors can be 
written in the form 

M? = WTt)T W) ( " 2ri + + 2S ' + 2 °> + r ' " 2S3) ' 

M » = 2(JTor(^ ( - 2A + A + M4 + M2 + A - 2 « • 

= 2(^)T(4 (2ri+, ' 4 " 2S4 " 2S2 + r3 + 2S3) ' 

M - = K^fr W) < 2ft + * " 2h ~ ^ + ft + 2 « • 

M 9 6 = 2 ^ + t)T (2ax + a 4 + 2 74 + 2 72 - « 3 - 2 74 ) , 
iN 

M W = 2(7T7)T (^ + 2 «3)- (83) 
3.3 Contribution from M 3 3 + M§ h 

These contributions need a treatment to some extend different from the preceding 
two subsections. First of all, we remove the factor Z introduced into M 3 ' 3 + 
M gh in Eq. (l53l) by means of (1381) and we introduce a separate notation for the 
corresponding contribution to the polarization tensor, 

nf A t sh = £ / dsdt e T (/ 2 , h 2 ) M 3 x 3 x ^ h (84) 

with 



M A 3 A Ui< = ^((((l-^) P -«)((l-^) P --))-trE^ 

A \ (A \ iF 



+E xy ^p + u)\^ p + u)+ tr-^ J . (85) 
Now we use the statements that the relations 

e T (/ 2 , h 2 ) (((1 - i) P - «) ((1 - P - «)) = - A e T (/ 2 , /> 2 ) (86) 

and 

Ot^ 2 , /> 2 ) ((^P + u) P + fi) + tr^ ) = -| 9 T (^, Z. 2 ) (87) 
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hold which can be verified by direct differentiation. In this way, in (1841) one 
integration can be carried out and we arrive at 



roo t—oo poo 

n 3 AA t gh = -E/ ds E xy e T (i 2 ,h 2 ) -E/ dts xy e T (i 2 ,h 2 " 

N Jo t =o N Jo 



s=0 



3.4 Tadpole and surface contributions 

In this subsection we collect the contributions resulting from the tadpole graphs 
given by Eq.( lT2~]) and the surface contributions which appeared in the preceding 
subsections. 

The tadpole contributions can be calculated easily since they are special cases 
of the basic loop contribution for s = collapsing the line of the charged gluon and 
keeping the line of the neutral gluon and for t = which the lines interchanged. 
All other rules remain valid so that these contributions can be written down 
easily, 

n tad P oi = sr[dt Q T (l 2 , h 2 ) ( 8\\> + A5 XX ) 

+ E / ds & T (l\ h 2 ) (ixE 8\ X > - 4sinh(2s)iF - (b\ x , + iFsinh(2s) + 5 
Now we collect the surface terms. A part of them has the form flBTl) . 



POD 

M surface = y dx Q n2 + 2 iF)M? x , 



(90) 



u=0 



with contributions to Mf x , from Eqs.([62]), fl65|). (j75|) and (j80l), 



M& = ~ (Pxx> - trE - R T XX , + 2QJ 1AA , - 2Qj 2xx , . (91) 

which can be rewritten in the form 



TTSurface 
ii AA' 



E I ds B T (l 2 , h 2 )M s n - £ [dt 9 T (/ 2 , h 2 )M s{ \ . (92) 

^ r J \t=0 AT J \s=0 



N " T - u N 



Now it can be checked that when adding (189]) and fl92|) and doing the obvious 
simplification all contributions cancel except for that which are proportional to 
u x u X i, i.e., to T^y, Eq.flSj). These collect into I1 D defined in Eq.ffTUl). 



n° = E / ds ^ 6 T (/ 2 , h 2 ) (-4 + 2tr£), -Ysfdt'Z ©t(/ 2 , h 2 ) (4 - 2tr£) 

N J Zs U=o N J Zt 



s=0 



(93) 
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Now we need from (12T]) and (J35 
6 T | 



6 



expj-jg) 
s exp 



Tl 



TV 2 \ 



isT 2 J 



sinh(s) (47r) 2 s 2 



iN 



(94) 



and after renaming the integration variables we obtain for the remaining contri- 
butions 



n 



D 



^ r°° dX iX 



N 



o A 2Tp4 



-46 



+ 46, 



(4tt) 2 



. , r 00 dA iN 
A 3 " 2Tp 4 



Acosh(2A) 
sinh(A) 



exp % 



N 2 



exp 



(95) 
1 



4: XT 2 J ' 



Thus, in our representation, the only contribution coming from the surface terms 
and tadpoles is the form factor U D . 



3.5 Debye mass of charged gluons 

The expression for I1 D derived in the preceding subsection can be a bit simplified 
by writing as a sum over N > (note the contribution from N = vanishes), 



IT 



(47T) 



E 

N=l 



dX sin(iV>4/T) Acosh(2A) 
A 3 " Tp 4 /N sinh(A) 



exp 



N 2 



4AT 2 J 



(96) 



It is obvious that for external momenta obeying p 4 = 2nlT (I integer) only the 
contribution from I = is nonzero and we arrive at the Debye mass in the field 
presence I1 D = —Si^m 2 ^, 



ml{B) 



1 

A^ 2 



E 

7V=1 



00 dX 
A 2 



cosh(2A) 
sinh(A) 



exp 



N 2 \ 

'axt 2 



(97) 



This expression coincides with the Debye mass of the neutral gluon, Eq.(123) in 
[3T] that may serve as a check of the correctness of carried out calculations. 

The integral over A is divergent at the upper limit due to the tachyonic mode. 
In fact, in this contribution one has to written the parametric representation from 
the very beginning on an axis rotated by 90 degrees clockwise and in the now 
divergent contributions rotate in anti-Wick direction as this was done in Ref 
As a result, in the high temperature limit, B/T 2 « 1, one obtains 



ml(B) 



-T 



-i 0.4775 




0.4775 



B 2 



(l6T 4 



(9£ 
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where the numeric values of the coefficients are substituted. We see that the real 
part in the magnetic background is smaller than without the field. The imaginary 
part is small because it appears in the next-to-leading order ~ yBT. 



4 Conclusion 

Here we collect the contributions which were calculated in the preceding subsec- 
tions. The form factors appearing in the decomposition ffTUl) of the polarization 
tensor read 

ITU 2 , h 2 + 2iF) = V / ds dt 9 T (/ 2 , h 2 + 2iF) Mi (99) 
n Jo Jo 



with 
Mi 



4(s + t) 2 -2(s-t) 2 cosh(2s) 
(I+t) 2 



M 2 = — \-2t 2 + A sinh(2s)t + 2 coslr (2s) + cosh(2s) [At 2 + 2 sinh(2s)t - 3j + 1 

M 3 = 1 (2(s + 5t) cosh 2 (2s) + At (t 2 + hst - 2) cosh(2s) 
4 f s I t ) _A 

+2(s + t) (2t 2 + 8 sinh(2s)t - 1 



M 4 = - 1 - - (A(s - t)t cosh 2 (2s) + (A(s - t)t + (s - 7t) sinh(2s)) cosh(2s) 
2(s + t)A v 



+8t(t - s) + (-2t d + Ust 2 + 7t- s) sinh(2s 
M 5 = ^-((cosh(2s)-l)(2t 2 + cosh(2s)-l)), 

M 6 = ^-((-2t 2 + cosh(2s)-l)sinh(2s)), (100) 

with A given by Eq. fl24|) . These are the same expressions as in our paper [29] 
dealing with the zero temperature case. 

The new temperature induced contributions from adding flTTj) and fl83|) read 

M iN A(s - 1) cosh(2s) 

7 ~ 2{s + t)Tp 4 7+t ' 

iN 4cosh 2 (2s) - 4 (2t 2 + 1) cosh(2s) 
M§ " 2(s + t)T P4 2A ' (101) 

iN tcosh 2 (2s) + (t + sinh(2s))cosh(2s)-2t + (2t 2 -l)sinh(2s)) 

Mq — — ; — — 4 : . 

9 2(s + t)T A 

Furthermore it holds M Wa = M 10 b = and Mu^ = — Mq. Here a comment is in 
order. The vanishing of Mio a and Mio& is a result of the calculations done here in 
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one loop order. At the moment it is not known whether there is some symmetry 
behind and whether this persists in higher loops. Similar remarks apply to the 
relation between Mg and M^y As a result, these two form factors contribute 
proportional to the tensor structure 

Tlx - All = u x id x , - id x u x > + iF(up) - ^y. (102) 

In addition to (1991) with the form factors (11001) and (11 Oil) we also have the 
contribution from the Debye mass f[9"B"j) which with ffTUl) and (jSJ) can be written 
in the form, 

U u„u v . (103) 

It should be mentioned that all the contributions (1991) . and, in particular, this 
one, are valid off the shell of the Matsubara values for the external momentum, 

1. e., for arbitrary j>4. In this case, of course, even the weak transversality does not 
hold and Eq. ( 1 1 3 f) is just the contribution on which this is realized. For p^ = 2nTl 
with integer I, from Eq.( l96l) I1 D = — Sifiin^B) follows and, as discussed in section 

2, the transversality holds. 

Finally we note that, like in the case of the neutral polarization tensor in [31], 
the form factors of the charged one contain imaginary parts. These result from 
the tachyonic mode and must be treated in the same way as in [31]. An example 
is Eq. (USD- 

One obvious application of the results obtained in our investigations is the 
resummation of perturbation series in the field at high temperature. That can 
be done by means of the solution of the Schwinger-Dyson equations for two-point 
Green's functions, where the derived tensor structures with arbitrary form factors 
have to be used as the input propagators. Whether or not the spectrum of gluons 
becomes stable due to the gluon magnetic mass, which is generated in the field 
at high temperature in this so-called super-daisy resummation, is an interesting 
problem for future investigation. 
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Appendix 

In this appendix we collect formulas which are necessary to restore the form 
factors from the expressions we are getting in the main text for the parts of the 
polarization tensor under the signs of the parametric integrals. In order not to 
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overload the notations we start from the T = case. A typical expression which 
appears from the calculation of the polarization tensor in the subsections 2.1 and 
2.2 has the form 

M = (Xp) x (Yp) x , + (Y T p) x (X r p) xi + W\\', (104) 

with 

X = r 5 11 + a iF + (3 5 ± , 
Y = s 5 11 +7 iF + 5 5 ± J 

W = a 5 11 +ciF + b 5 ± . (105) 

We note that the parameters in fll05j) has to fulfill certain relations if M shall 
obey the weak transversality condition ([1]). Irrespective of that, the first four 
form factors can be restored. Using (jl]) we rewrite (11041) in the form, 

M = -2rsT (1) -2(a 7 + /35)T (2) - (rd + s(3)T {3) + (r 7 - sa)T (4) 

+ A5 ll +CtF + B5 ± (106) 

with 

A = 2rsl 2 + (r8 + s(3)h 2 + r 7 - sa + a , 
B = 2(35h 2 + (r5 + s(3)l 2 + ad - /3j + b , 

C = (sa - n)l 2 + (a5 - (3-f)h 2 + 2a 7 + 4(35 + c . (107) 

From (11061) . the first four form factors read 

M 1 = -2rs, M 2 = -2{aj + P6), 

M 3 = -(rd + sfi), M 4 = n-sa. 1 J 

In the main text, the parts of the polarization tensor appear (partially after 
integrating by parts) in specific forms matching one of the following representa- 
tions of W, 

W x = -2Y x \,(pXp), 
W 2 = -2X X y(pY^p) , 

W 3 = -(Wi + W 2 ). (109) 

First we consider W\ and use (pXp) = rl 2 — a + f3(h 2 + 2iF) to obtain 

W l = -2 (s 5 11 - 7 iF + 5 5 X ) (rl 2 -a + (3(h 2 + 2iF)) , 

= -2s ( (rl 2 -a + f3h 2 ) 5 n + (-7 (rl 2 -a + (3h 2 ) + (35) iF 

+ (5(rl 2 -a + (3h 2 )-(3^)5 ± ), (110) 
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from which we identify the corresponding expressions for a, b and c in the second 
line in Eq. (ll05j) . These we insert into A, B and C in (11071) . 

A = (r5 — sf3)h 2 + + sa , 

B = -(r5 - sf3)l 2 + 3(a5 + (3^) , 

C = (r 7 + sa)l 2 + (ad + P^)h 2 . (Ill) 

Comparison with (j3J) shows that these contributions collect just into form factors 
such that M, (11041) . with W\ from (1 11 01) is weak transversal. From fllllD and (jlj) 
the remaining form factors follow, 

M 5 = r5-sp } M 6a = r 7 + sa, M 6b = ad + fa. (112) 

The same procedure can be repeated for W%. The result is just that the form 
factors M 5 , M 6a and M 6 b are the same but with reversed sign. Finally, for W3, 
they just compensate each other. Collected together, these results read 





M 5 


M 6a 


M 6b 


Wi 


r5 — s/3 


r 7 + sa 


a5 + /?7 


w 2 


-(r6-sP) 


— (rj + sa) 


-{aS + fii) 


w 3 












(113) 



Also we need the special case when Y T = X. In that case there is only one 
expression for W and we have to consider (this is M from (1104p divided by two) 

M=(Xp) x (x T p) y + W xy , (114) 

with 

W = -X xx ,(pXp) , (115) 

In this case the nonzero form factors are 

Mi = -r 2 , M 2 = al 2 - p 2 , M 3 = -r/3, M 4 = -rat . (116) 

Now we generalize these formulas to finite temperature. The basic expressions 
which appear in the subsections above have the form 

M=(Xp- fiu) x (Yp - vu) x , + (Y T p - vu) x (X T P - fiu) y + W xx >, (117) 

Here \i and v are numbers and for W again 3 types of expressions appear, 

Wi = -2 (Vv " ((P^P) " A*(«P)) , 

w 2 = -2 (^x AV - /i^) ((pi^p) - 1/ w) , 

W 3 = \(W 1 + W 2 ). (118) 
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Let us consider the contributions AW which are new as compared to Eg. fl!04|) . 
We start from the easiest part, namely that proportional to \xv. As easily can be 
seen from ( 11171) and ( 11181) they cancel so that only contributions linear in \i and 
linear in v remain. For W\ these read 



AWi 



—v 



-II 



(up) 



• s r T (7) _ rp{10a)} _j_ . g . ^(8) _ r (106)^ _j_ 7T (9) 



(119) 



(ixp) ^ ' (up) 

Here we used in addition to (jjTj) the following formulas which are easy to verify. 



and 



^(7) + y(10a) 
T (7) _ T (10a) 
T (8) + T (10b) 
rp(T) _ rp(10b) 



rp(9) _ 2J"(*) 



(u a /a' + /a«A') (up) - 2u x u\d 2 , 
(m a /a' + huy) (up) - 25 a ' a ,(mp) 2 , 
(u\h\> + h x u y ) (up) - 2u x u y h 2 , 
(w a /ia' + h x uy) (up) - 25 xx ,(up) 2 



u x idy - id\U\> 



(up) 



U\U X > 



(120) 



:i2i) 



with 



In a similar way we get 



AWo 



—v 



T« = _ ( UAtiA , + («p) 2 iF) 

= t^t ~~ o — ~ f (upfT^ - T^) 

(up) p 2 + 2iF v v yj J 

P 



(122) 



(up) 



(up) \ ' (-up) 

Taking the half sum of both the expression 



( T (7) _ T (10a)) + ^(8) _ T (106)^ _ aT (9) (123) 

^(7) + ^(lOtt)^ + ° ^ T (8) + T (lOb)^ + ^ ^(9) _ 2T ( *^ 



aw* 



T (7) + ^L T (8) _ a f T 0) _ T (*: 
(?/.rn v 



(up) 



(up) 



-ft 



(up) 



y(7) + _|_ 7 M9) _ T (*: 

(up) \ 



(124) 



emerges. 
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Finally we note the special case with Y = X T , fi = v for which we have 

M=(Xp- fiu) x (x T p - fiu) v + Wxx>, (125) 



with 

Ufjfll V 



W = -( X\\i - ((pXp) - M H) . (126) 



(up) 

We obtain in addition to (jllfip 



AW = -11 



(ii.n\ V 



(*) 



(up) (up) 
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